Slow-roll k-essence 
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We derive slow-roll conditions for thawing k-essence with a separable Lagrangian p {X, (/>) = 
F{X)V{<j))- We examine the evolution of the equation of state parameter, w, as a function of 
the scale factor a, for the ceise where w is close to —1. We find two distinct cases, corresponding 
to X « and Fx ~ 0, respectively. For the case where X m the evolution of ((> and hence 
w is described by only two parameters, and w{a) is model-independent and coincides with similar 
behavior seen in thawing quintessence models. This result also extends to non-separable Lagrangians 
where X « 0. For the case Fx ~ 0, an expression is derived for u){a), but this expression depends 
on the potential V{<j)), so there is no model-independent limiting behavior. For the X « case, we 
derive observational constraints on the two parameters of the model, wq (the present-day value of 
w), and the K, which parametrizes the curvature of the potential. We find that the observations 
sharply constrain wo to be close to —1, but provide very poor constraints on K. 

PACS numbers: 98.80.Cq ; 95.36.-l-x 



I. INTRODUCTION 

Cosmological data from a wide range of sources includ- 
ing type la supernovae P, 0, Q , the cosmic microwave 
background Q, baryon acoustic oscillations @, @L dus- 
ter gas fractions 0, [1] and gamma ray bursts [3, seem 
to indicate that at least 70% of the energy density in the 
universe is in the form of an exotic, negative-pressure 
component, called dark energy. (See Ref. [ll[ for a re- 
cent review). 

The dark energy component is usefully parameterized 
by its equation of state (EoS) parameter w, defined as the 
ratio of its pressure to its density. Observations constrain 
w to be very close to —1. For example, if w is assumed 
to be constant, then -1.1 < w < - 0.9 

On the other hand, a variety of models have been pro- 
posed in which w is time varying. A common approach 
is to use a scalar field as the dark energy component. 
The class of models in which the scalar field is canonical 
is dubbed quintessence [13, [H, [H, and has been ex- 
tensively studied. In most of these models the field rolls 
slowly on a flat potential. It is also possible to have mod- 
els in which the field is trapped in a false vacuum (see 
[ist and references therein for examples of such models). 

A related, yet somewhat different approach is phantom 
dark energy, i.e., a cornponent for which w < —1, as first 
proposed by Caldwell Il9l . Such models have well-known 
problems [23, [2lj, [l^, l23j | (however see for recent at- 



tempt to construct a stable model) , but nevertheless have 
been widely studied as potential dark energy candidates. 

In both the above approaches, the potential energy of 
the scalar field is responsible for bringing about the late- 
time acceleration of the Universe. A third type of model 
utilizes the kinetic energy of the field through the use 
a non-canonical kinetic term in the Lagrangian. Such 
fields were first used in the context of inflation, in a sce- 
nario that is called k-inflation [2^ . They have since been 
studied as dark energy candidates [2^, 21, 2M and these 
models are called "k-essence" . 

Given the considerable freedom that exists in choos- 
ing the potential function of the scalar field V ((/>) , as 
well as the kinetic function F {X) for k-essence, it would 
be interesting to find any model-independent evolution 
for w{a). Some recent work has uncovered such model- 
independent evolution for quintessence and phantom 
models when the field evolves in the vicinity of extrema 
of the potential, and w ^ —1. Ref. [l^ considered the 
evolution of a a scalar field, initially at rest, in a potential 
satisfying the "slow-roll" conditions: 
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The first condition ensures that w is close to — 1, 
while the two conditions taken together indicate that 
{\/VMdy/d4)) is nearly constant. In the terminology of 
Ref. [301, these are "thawing" models. 

For all potentials satisfying these conditions, it was 
shown in [29| that the behavior of w can be accurately 
described by a unique expression depending only on the 
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present-day values of fi^ and the initial value of w. In 
[3ll | this result was extended to phantom models satisfy- 
ing Eqs. (IlEI), and the w dependence of these phantom 
models was shown to be described by the same expression 
as in the quintessence case. 

The slow roll conditions. Eqs. ((lEl); while sufficient 
to ensure w ~ — 1 today, are not necessary. In [s^l, a 
second possibility was considered, in which equation ([1]) 
holds, but equation ^ is relaxed. This corresponds to 
a quintessence field rolling near a local maximum of its 
potential. As in the case of slow-roll quintessence, this 
case can be solved analytically. In this case, there is an 
extra degree of freedom, the value of {l/V){(PV/d(t)'^), so 
that instead of a single solution for the evolution of w, one 
obtains a family of solutions that depend on the present- 
day values of fi^ and w and the value of {l/V){(PV/d(j)'^) 
at the maximum of the potential. This family of solutions 
includes the slow-roll solution as a special case in the 
limit where {l/V){cPV/d4>'^) 0. The corresponding 
result for phantom fields, where the field rolls near the 
minimum of its potential was derived in (ssj - and an 
identical expression for w{a) was obtained. In both the 
above cases, it was shown that the expression for w{a) 
reduces to the corresponding ones in [23] and [sH as the 
potential gets flat, i.e., {l/V){d'^V/d(l}'^) 0. 

Ref. [131 examined the opposite case of dark energy 
models in which a quintessence or a phantom field rolls 
near the vicinity of a local minimum or maximum, re- 
spectively, of its potential. It was shown that as long as 
Eq. Il]) is satisfied, (although Eq. ([2]) need not be), the 
evolution of w is described by an expression identical to 
the one in [s^, [13] . In these cases the evolution of w was 
found to encompass a richer set of behaviors, including 
oscillating solutions. 

In [1^, one of us (TC) showed that the expression de- 
rived in [13, [13, [13| has a wider applicability than the 
cases of fields rolling close to extrema in their potentials. 
The following more general slow-roll conditions on the 
potential are derived there: 



e < 1, 



6HW 
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(3) 
(4) 



and while retaining the assumption that w ~ —1, 
dropped the assumption that the field is close to a lo- 
cal extremum in the potential. Interestingly, once again 
the expression for the evolution of w under these more 
general conditions was found to coincide exactly with the 
one in refs. [H, [H, [13| . 

In this paper we investigate thawing models in k- 
essence. We derive slow-roll conditions for thawing k- 
essence analogous to equations ©-dH), and show that, 
when w ~ —1, in some cases one does obtain the same 
model independent evolution of w seen in the above ref- 
erences. 



II. SLOW-ROLL THAWING K-ESSENCE 

The Lagrangian density of k-essence [2^, [13, [13] is 
p{(f>,X), where X — — V^(/!)V^0/2. The pressure of 
the scalar field is given by X) itself and the energy 
density is given by — 2X{dp/dX) — p [H, [13], so 
that the equation of state parameter, w, is 



w 
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2X{dp/dX)~p 
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Working in units of 8ttG = 1, the basic equations in a 
flat universe are 
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((1 + 3wb)pb + (1 + 3w)p0) , 



(6) 



(7) 



(8) 



where pB and pB are the energy density and the pressure 
of the background matter and/or radiation, respectively. 

We now proceed to derive slow-roll conditions for k- 
essence with the following factorized form of AT): 



p{^,X) = V{^)F{X). 



(9) 



The equation of motion of the scalar field is then written 
as 

^ (Fx + 2XFxx) + iHFx^ + {2XFx - F) ^ = 0, 

(10) 

where Fx = dF/dX and V = dV/d(j). We also intro- 
duce the sound speed of k-essencc, which is the relevant 
quantity for the growth of density perturbations. 



2 _ dp/dX 



Fx 



dp/dX 2XFxx + Fx ' 



(11) 



Using Cg, The equation of motion Eq. pUj) is rewritten 
as 



,2XFx-FV' 



Fx 



V 



= 0. 



(12) 



A. Slow-Roll Conditions for K-Essence 

By slow-roll k-essence, we mean a model of k-essence 
whose equation of state w is close to —1 so that 



\XFx\ « \F\ 



(13) 
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Thawing models correspond to the equation of state w = 
Pcjil P<j> very close to —1, so that the Hubble friction is not 
effective and hence is not necessarily small compared 
with 3iJ(/) in Eq. ((1 
We derive the 
k-essence during 
epoch. Generalizin; 
quintessence [ssl. Iscl. 
function: 



slow-roll conditions for thawing 
the matter/radiation dominated 
j; the corresponding expression for 
|37| . we first introduce the following 



(14) 



As stated above, for thawing models, /? is a quantity of 
0(1). We assume f3 is approximately constant in the 
sense that \$\ ^ H\(3\, and the consistency of this as- 
sumption will be checked later. In terms of (3, from Eq. 
(fT2l) using Eq. (fT3|) . (p is written as 



FV 



3{l + f3)HFxV' 



and the slow-roll condition Eq. p3)) becomes 



\F\V'- 



6m\Fx\V^ 



«1, 



(15) 



(16) 



where we have omitted 1 + /3 since it is an 0(1) quan- 
tity, and we have introduced the factor of 1/6 so that 
e coincides with the slow-roll parameter for thawing 
quintessence [H], e ^(V'^/H^V). ^ Eq. ([HI) is a 
k-essence counterpart of the quintessence slow-roll con- 
dition V'^/H^V < 1. 

Similar to the case of inflation, the consistency of Eq. 
(fT4)) and Eq. p2|) should give the second slow-roll con- 
dition. In fact, from the time derivative of Eq. p5|) we 
obtain 



'V~~V^) 9(1 + P)H^Fx 



(I + wb) 



(17) 



where we have used H / ~ — 3(1 + 1(75) /2 from Eqs. 
d?])-® and have assumed (3 < Hj3. While the left-hand- 
side of Eq. p7|) is an almost time-independent quantity 
by assumption, the first term in the right-hand-side is a 
time-dependent quantity in general. Therefore the equal- 
ity holds if the first term is negligible, which requires in 
addition to Eq. (HI 



FV" 



iH^FxV 



so that (3 becomes 



a- 



(18) 



(19) 



or rj itself becomes a constant so that 

^ = -3(l + /3) (p-\{1 + wb) 



(20) 



The former condition would correspond to the slow-roll 
models with X ~ 0, while the latter corresponds to the 
slow- roll models with Fx — 0. The expression for (3 
given by Eq. (|19p is approximately constant, which is 
consistent with our assumption. Here the factor —1/3 
is introduced in Eq. ((T8|) so that 77 coincides with the 
slow- roll parameter for thawing quintessence (35| . 77 = 
^{V" / H"^). ^ Eq. (fT5)) is a k-essence counterpart of the 
quintessence slow- roll condition \V"\/H'^ <C 1. 

Eq. ([TB)) and Eq. (|18p constitute the slow-roll condi- 
tions for thawing k-essence during the matter/radiation 
epoch. Note that these expressions have assumed a neg- 
ligible contribution to the expansion rate from the k- 
essence itself, and so /3 is no longer a constant and Eq. 
([19]) or Eq. ([20|l becomes progressively less accurate 
as the k-essence begins to dominate at late times. In 
what follows, we do not make the assumption of mat- 
ter/radiation domination, so that our results will be ac- 
curate up to the present. 



B. Parametrizing the Equation of State 

Next we derive general solutions for in the limit 
where |1 -|- wj ^ 1, and we derive w as a function of 
a. We note that 1 -|- w = implies {&) X = Q or (b) 
Fx = 0. In the following we consider each case. 

Case (a): First we consider the case where X ~ 0. In 
this case, w 1, and Eq. (fT^ simplifies to 



^ Fx{0)V 



(21) 



The Hubble friction term in Eq. (I12[) can be eliminated 
by the following change of variable [34I 



.3/2 



(22) 



where (pi is an arbitrary constant, which is introduced for 
later use, and then Eq. (fT2|l becomes 



u+ -;{PB +P4,)i 



3/2 ^(0)^' 
Fx{Q)V 



0. 



(23) 



We assume a universe consisting of matter and k-essence 
with ui ~ — 1. Then the pressure is well approximated 
by a constant: ps + p<p ~ p<j) — —P4,o, where p^o is the 



^ In terms of a canonically normalized scalar field ip, 3r]H^ 
1 A canonically normalized scalar field ip corresponds to dtp = + 5^?=/^. and when combined with e < 1, r; 

yl/2rf0 for F = X - 1. i{V,y,y,/H'^). 
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nearly constant density contributed by the k-essence in 
the hmit w ~ — 1. Eq. then becomes 



u - -p4,ou + a' 



3/2 



F{0)V' 
Fx{0)V 



= 0. 



(24) 



Since we consider a slow- roll scalar field {X ~ 0), the 
potential may be generally expanded around some value 
which we identify with the initial value, in the form 
(up to quadratic order) 



F(0) = V^((/),;) + ^^'(^O(0-0.) + 



Substituting the expansion given by Eq. 
((M)) and taking p^o = -FiO)V{(t>t) gives 



FiO)V"{4'^) 
FxiO)V{(b^) 



-FiO)V{cj,,)ju^ 

Fx(o)v(4>,) : 



(25) 
into Eq. 

(26) 



This equation is identical to Eq. (19) in Ref. [3a| by 
the substitution: V -FV,V" -> ~FV" /FxV.V 
~FV' / FxV . Therefore, the evolution of is the same 
(in functional form) and the equation of state i s ag ain 
given by the same functional form derived in [H, l35j |: 



1 + w{a) = (1 + wo)a' 



3(A'-i) ( (K F{a)){F{a) + 1)^-^ + {K + F{a)){F(a) 1)^' 



(27) 



where K and F{a) [not to be confused with F{X)] are 
defined by 



K 




iFx{Q)V{<Pi 



F{a) = J I + (r!7„i - l)a 



(28) 
(29) 



This equivalence with a scalar field having a canonical 
kinetic term can also be easily seen by noting that F{X) 
can be expanded (if it is analytic at X = 0) for small X 
as F[X) = i^(0) + Fx(0)X and the Lagrangian reduces 
to that of a canonical scalar field by field redefinition. 
This implies that the equivalence is not limited to the 
factorizcd form of p{(j),X) (Eq. [9]), but holds for more 
general p(0, X) because p{(f>, X) can be expanded (if it is 
analytic at X = 0) for small X as p{(j), X) ~ p{(j), 0) + 
Px{cl),0)X. 

Case (h): For the case of Fx ~ 0, we may expand 
F{X) around the extremum of F (say aX X = Xm) as 
(like hilltop quintessence) 



F{X) = F{X„ 



-^Fxx{Xm){X — Xm) 



(30) 



In this case, assuming X ~ X^, Eq. simplifies to 



X + 3H{X - Xm) 



F{Xm) 



d 



2X„iFxx{Xm) dt 



{\nV) = 0. (31) 



Therefore, the evolution of X depends on the shape of 
V{(l)), which is no longer Taylor expanded because 4> can 
evolve significantly in this case. 

However, one can still derive an analytic expression for 
w{X) . Starting from Eq. (O, one can write w{X) as 



w{X) 



F{X) 



2XFx {X) - F{X) ■ 



(32) 



For Fx{X) 0, this expression can be expanded about 
the point X = Xm- In this case terms up to second order 
must be retained, because while {X — X,„)^ is small, the 
second Taylor coefficient can be large leading to a large 
contribution. The resulting expansion gives: 



w{X) = -l- 



2X„iFxx{Xm) 

F{Xm) 



{X - Xm) - 



2Fxx{Xm) 

F{Xm) 



^^m^Xxi^rn) 
F^Xm) 



(X — XmY 



(33) 
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Our numerical results indicate that equation ([55)1 is an 
excellent approximation for the evolution of w for the 
case where Fx ~^ (see Fig. 4 below). However, it is 
of limited usefulness, since 'w{a) in this case depends on 
X{a), and X{a), in turn, depends on the functional form 
of V^(0). 

III. COMPARISON TO NUMERICAL RESULTS 

We now turn to numerically solving the equations of 
motion in order to compare Eq. (j27p against the exact 
evolution. We consider a Universe consisting of perfect 
fluid dark matter and k-essence dynamical dark energy. 

We work with models which satisfy the slow roll condi- 
tions derived above, i.e., Eq. and Eq. For the 
cases where X « 0, the three specific models considered 
are listed below. 

1. Case 1: 

F{X) = VI - mX, V (0) = ^e-'^'/'"', see Fig. (P. 
This is the rolling tachyon Lagrangian [38|1 sug- 
gested by the boundary string field theory [sgj . 

2. Case 2: 

F{X) = VI - mX, V{(j)) = A(j)-°', see Fig. ©. 
This is the model studied in [40l |. 

3. Case 3: 

p(<^,X) = mT(^)-mT(^)(l-2X/T(</)))i/2_F(0), 
T{(j)) = (j)*, V{(f>) = ix^cj)"^ see Fig. ©. This is the 
Dirac-Born-Infeld (DBI) model discussed in Refs. 

ESI 

where m = ±1 for all three cases. In the first two cases, 
we choose 4>{t = ti) ~ 1 and (f){t = ti) = 0. The constants 
in the potentials {A, a and a) are then adjusted to give 
P<t> ~ Pa ^-nd w = —1 at t = ti and w = —0.9 or w = —1.1 
at t = to- For the third case, we choose 4>{t = ti) = 
and adjust cpi and /x to get the above initial and final 
conditions. 

In all these cases we find excellent agreement between 
the numerical and analytic results, i.e., Sw/w < 0.01. 
The success of our approximation for the DBI case in- 
dicates that, as noted earlier, it is not just restricted to 
separable Lagrangians. Note that for the second case it 
is found in Ref.[40[ that models have a unique equation 
of state, which corresponds to = 1 in Eq. (|27|1 . 

As noted earlier, for the set of models where Fx ~ 0, 
we do not expect the evolution of w{a) to converge to a 
common potential-independent behavior. This is evident 
from Fig. (j4]) where we plot the 'w{a) behavior for two dif- 
ferent potentials, V {(jj) = Ae"'^'/'^' and V {(j}) = Af/)"", 
and in both cases we choose F (X) = X^ + {X — Xj^f ■ 
We take initial conditions 0, = 1 and X ~ Xm and 
adjust the constants appropriately. The w{a) behaviors 
turn out to be very different, indicating that there is 
no common model-independent behavior for this class of 
models. On the other hand, our analytic approximation 
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0.2 0.4 0.6 0.8 1 

a 

FIG. 1: Evolution of ■w{a) for a k-essence model with 

F{X) = VI -mX and V {(j)) = Ae"*''' ' "'^ . The upper plot 
corresponds to tlie m = 1 case and the lower one corresponds 
to the m = —\ case. The solid (blue) lines denote the nu- 
merical result and the broken (red) lines denote the analytic 
approximation given by Eq. ([27]). 

-0.85 1 ' ' ' ' 1 




0.2 0.4 0.6 0.8 1 

a 

FIG. 2: Evolution of wia) for a k-essence model with F (X) = 
Vl — mX, V (0) = The upper plot corresponds to the 

m = \ case and the lower one corresponds to the m = — \ case. 
The solid (blue) lines denote the numerical result and the 
broken (red) lines denote the analytic approximation given 
by Eq. (t27)l . 



(equation [33]) does give excellent agreement, although of 
course the approximation is itself a function of ^(0) in 
this case. 



IV. OBSERVATIONAL CONSTRAINTS IN THE 
(it)o,A') PLANE 

The results of this paper, combined with previous stud- 
ies, indicate that Eq. ([?7)) applies both to quintessence 
models and to a subset of k-essence models with ~ — 1. 
Hence Eq. ([27]) is a useful and physically well-motivated 
parametrization for w[a) that can be compared with the 
observations. So, in this section, we present the obser- 
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FIG. 3: Evolution of w{a) for a DBI model with p (0, X) = 

mr(0)-mr(0)(i-2X/r(0))i/2_v/(0), r(</<) = </>^ 1^(0) = 

fj.^(f>^. The upper plot corresponds to the m = 1 case and the 
lower one corresponds to the m = — 1 case. The solid (blue) 
lines denote the numerical result and the broken (red) lines 
denote the analytic approximation given by Eq. (|27p . 
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w_0 



FIG. 5: Contours at 68.3% (red, inner), 95.4% (blue, mid- 
dle), 99.7% (green, outer) confidence level on wq and K. The 
Constitution SN set was combined with BAO constraint. 
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"0 0.2 0.4 0.6 0.8 1 



a 

FIG. 4: Solid lines give evolution of 'w{a) for the models in- 
dicated, for both of which Fx ~ throughout the evolution. 
Broken lines give the analytic approximation. Both the ex- 
act evolution and the analytic approximation depend on the 
choice of V{(j>), as shown. 



lations (BAO) from the SDSS data [1]. ^ Uncertainties 
in the distance modulus of a supernova include uncer- 
tainties in light curve fitting parameters (the maximum 
magnitude, stretch parameter, color correction parame- 
ter )and due to the peculiar velocity (400 kms^^) as given 
in g. 

BAO measurements from the SDSS data provide a con- 
straint on the distance parameter A defined by 

M'>-(n..Biy"{jj^,£^)"' (34, 

to be A{z = 0.35) = 0.469 ± 0.017.'^ 

The joint constraints from SNIa and BAO are shown 
in Fig. m We marginalize over flm to derive the con- 
straints. The allowed range of wq is narrow: —1.04 < 
Wo < — 0.86(1(t). ^ We find that the cosmological con- 
stant Wo = —1 is fully consistent with the current data. 
Note that K, which parametrizes the curvature of V{(t)i), 
is not well-constrained by current SNIa and BAO data. 



vational constraints on the equation of state parameters 
Wq and on K . 

First, we note that the cosmological constant corre- 
sponds to a line in the {wo,K) plane: wq = —1 irrespec- 
tive of K . This can be understood for a canonical scalar 
field by noting that wq — ~1 corresponds to the case 
where the scalar field sits at the minimum (K < 1) or 
the maximum {K > 1) of the potential. 

As observational data we consider the recent compila- 
tion of 397 Type la supernovae (SNIa), called the Con- 
stitution set with the light curve fitter SALT, by Hicken 
et al. 0] and the measurements of baryon acoustic oscil- 



V. CONCLUSIONS 

Our results indicate that k-essence models with w near 
— 1 can be divided into two broad categories: models with 



^ We do not consider BAO distance measurements by Percival et 
al. Q because some points of tension were noted between the 
data sets Qj. 

Note that A changes only slightly in varying the spectral index 
of the matter power spectrum, n. This value is for n = 0.98, 
which changes to A = 0.472 if n = 0.96. 
^ We note that SNIa data alone do not constrain uiq much; —1.2 < 
wo < -0.7. 
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X « and those with Fx ~ 0. In the former case, we 
find a generic evolution for w{a) which is identical to 
the previously-derived evolution for quintessence. This 
strengthens the case that equation (|27p is a useful and 
physically well-motivated parametrization for u;(a) that 
can be compared with the observations, since it applies 
both to quintessence models and to a subset of k-cssence 
models with w k. —1. Applying this parametrization to 
SNIa data and BAO, we find that the present-day value 
of w is constrained to lie near —1, while the curvature 
parameter K is poorly constrained by the observations. 
Further, we see that the cosmological constant limit of 
these models is consistent with the current data. 

On the other hand, k-essence models with Fx ~ can 
demonstrate quite different behavior. In this case, the 



evolution of w{a) is strongly dependent on the particular 
potential, and there is no "generic" behavior. 
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